Nonlinear transport and oscillating magnetoresistance in double 

quantum wells 

N. C. Mamani\ G. M. Gusev\ O. E. Raichev^-^, T. E. Lamas\ and A. K. Bakarov^'^ 
^Instituto de Fisica da Universidade de Sao Paulo, 
CP 66318 CEP 05315-970, Sao Paulo, SP, Brazil 
"^Institute of Semiconductor Physics, NAS of Ukraine, 
Prospekt Nauki 45, 03028, Kiev, Ukraine 
^Institute of Semiconductor Physics, 
Russian Academy of Sciences, Novosibirsk 630090, Russia 
(Dated: November 27, 2008) 

Abstract 

We study the evolution of low-temperature magnetoresistance in double quantum wells in the 
region below 1 Tesla as the applied current density increases. A flip of the magneto- intersubband 
oscillation peaks, which occurs as a result of the current-induced inversion of the quantum com- 
ponent of resistivity, is observed. We also see splitting of these peaks as another manifestation 
of nonlinear behavior, specific for the two-subband electron systems. The experimental results 
are quantitatively explained by the theory based on the kinetic equation for the isotropic non- 
equilibrium part of electron distribution function. The inelastic scattering time is determined from 
the dependence of the inversion magnetic field on the current. 

PACS numbers: 73.23.-b, 73.43. Qt, 73.50.Fq 
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I. INTRODUCTION 



The nonlinear transport of electrons in two-dimensional (2D) electron systems placed 
in a perpendicular magnetic field has been extensively studied in the past in connection 
with the breakdown of the quantum Hall effect at high current densities.^ More recently, it 
was realized that the current causes substantial modifications of the resistance even in the 
region of weak magnetic fields and relatively high temperatures, when the Landau levels are 
thermally mixed so the Shubnikov-de Haas oscillations (SdHO) are suppressed. 

The present interest to the static (dc) nonlinear transport in 2D systems is stimulated by 
observation of two important phenomena. First, in high-mobility systems there appears a 
special kind of magnetotransport oscillations, when the resistance oscillates as a function of 
either magnetic field or electric current. Second, it is found that the current substantially 
decreases the resistance even at moderate applied voltages. ^'^ The observed phenomena are of 
quantum origin, they are caused by the Landau quantization of electron states and refiect the 
infiuence of the current on the quantum contribution to resistivity. The oscillating behavior 
is explained by modification of the electron spectrum in the presence of high Hall field, ^'^'^ 
while the decrease of the resistance is most possibly governed by modification of electron 
diffusion in the energy space, which leads to the oscillating non-equilibrium contribution 
to the distribution function of electrons.^ A theory describing both these phenomena in a 
unified way has been recently presented.^ 

In contrast to the Hall field-induced resistance oscillations, the phenomenon of decreasing 
resistance has not been studied extensively in experiment. Though the available data^ 
support the theory^'^ predicting nontrivial changes in the distribution function as a result 
of dc excitation under magnetic fields, they are not sufficient for definite interpretation of 
the observed phenomenon in terms of this theory. For better understanding of the physical 
mechanisms of nonlinear behavior, further investigations are necessary. 

In this paper, we undertake the studies of nonlinear magnetotransport in double quantum 
wells (DQWs), which are representative for the systems with two closely separated occu- 
pied 2D subbands. In contrast to the quantum wells with a single occupied subband, the 
positive magnetoresistance,^ which originates from the Landau quantization, is modulated 
in DQWs by the magneto-intcrsubband (MIS) oscillations.^'^ These oscillations, whose max- 
ima correspond to integer ratios of the subband splitting energy A12 to the cyclotron energy 
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hwc^ are caused by periodic variation of the probability of elastic intersubband scattering 
of electrons by the magnetic field as the density of electron states becomes an oscillating 
function of energy. As a result, the changes in the quantum contribution to the conductiv- 
ity are directly seen from the corresponding changes of the MIS oscillation amplitudes. In 
particular, we observe a remarkable manifestation of nonlinearity in DQWs, the inversion 
of the MIS oscillation picture, which appears when the quantum magnetoresistance changes 
from positive to negative as a result of increased current (Fig. 1). By adopting the ideas of 
the theory of Ref. 7, we explain basic features of our experimental data and determine the 
inelastic relaxation time of electrons in our samples. 

The paper is organized as follows. In Sec. II wc describe the experimental details and 
present the results of our measurements. In Sec. Ill we generalize the theory of Ref. 7 to 
the case of two-subband occupation. A discussion, including comparison of experimental 
results with the results of our calculations, is given in Sec. IV. The last section contains the 
concluding remarks. 

II. EXPERIMENT 

The samples arc symmetrically doped GaAs double quantum wells with equal widths 
dw = 14 nm separated by AlajGai-^^As barriers with width db=1.4, 2, and 3.1 nm. Both 
layers are shunted by ohmic contacts. Over a dozen specimens of both the Hall bars and van 
der Pauw geometries from three wafers have been studied. We have studied the dependence 
of the resistance of symmetric balanced GaAs DQWs on the magnetic field B at different 
applied voltages and temperatures. While similar results has been obtained in all samples 
with different configuration and barrier width, we focus on measurements performed on two 
samples with barrier width dh=lA nm. The samples have mobilities of 9.75 x 10"'' cm^/V s 
(sample A) and 4.0 x 10^ cm^/V s (sample B) and total electron density Ug — 1.01 x 10^^ 
cm~^. The samples are Hall bars of width 200 //m and length 500 //m between the voltage 
probes. The resistance R — R^x was measured by using the standard low-frequency lock-in 
technique for low value of the current. We also use DC current, especially for high-current 
measurements. The results obtained with AC and DC techniques are similar. The subband 
separation A12, found from the MIS oscillation periodicity at low B, is 3.7 meV for sample 
A and 5.1 meV for sample B. 
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FIG. 1: (Color online) Magnetoresistance of the sample A for three different currents / at T = 1.4 
K. The oscillations are inverted with the increase of the current. The inset shows the linear and 
non-linear (at / = 200 /uA) magnetoresistance in the low-field region. 

The resistance of the sample A as a function of magnetic field at different currents is pre- 
sented in Figs. 1 and 2. At small currents, the magnetoresistance is positive and modulated 
by the large-period MIS oscillations clearly visible above B = 0.1 T. The small-period SdHO, 
superimposed on the MIS oscillation pattern, appear at higher fields in the low-temperature 
measurements (Fig. 1). With increasing current /, the amplitudes of the MIS oscillations 
decrease, until a flip of the MIS oscillation picture occurs. This flip, which we associate with 
inversion of the quantum component of the magnetoresistance from positive to negative, 
starts from the region of lower fields and extends to higher fields as the current increases. 
Therefore, one can introduce a characteristic, current-dependent inversion field Binv The 
inset to Fig. 2 shows the behavior of the magnetoresistance near the point of inversion. In 
this point, apart from the transition from negative to positive quantum magnetoresistance, 
we observe an additional feature that looks like splitting of the MIS oscillation peaks or 
appearance of the next harmonic of the MIS oscillations. This feature persists in higher 
magnetic fields. In contrast to the MIS oscillations, the SdHO are not inverted by the cur- 
rent, as seen in Fig. 1. However, the SdHO amplitudes decrease as the current increases 
until the SdHO completely disappear in the low-field region. We attribute this suppression 
of the SdHO to electron heating at high current densities. 

The amplitudes of inverted MIS oscillations increase with increasing current and become 
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FIG. 2: (Color online) Magnetoresistance of the sample A for different currents at T = 4.2 K. The 
inset shows inversion of the quantum magnetoresistance around B = 0.2 T. 

larger than the MIS oscillation amplitudes in the linear regime. At low temperatures the 
ratio of the corresponding amplitudes varies between 2 and 3; see Fig. 1. However, when 
the current increases further, the amplitudes of inverted peaks slowly decrease, this decrease 
goes faster in the region of lower magnetic fields. This property is seen in Figs. 3 and 
4, where the magnetoresistance data for the sample B is presented. The typical current 
dependence of the inverted peak amplitudes at T = 1.4 K is shown in the inset to Fig. 3. 
The behavior of magnetoresistance at 4.2 K, shown in Fig. 4, is similar. In the chosen 
interval of magnetic fields, the SdHO at 4.2 K are suppressed even in the linear regime. The 
splitting of the MIS oscillation peaks is clearly visible in Fig. 4 at / = 80 /iA. For / = 100 /iA 
this splitting apparently develops in the frequency doubling of the MIS oscillations. Further 
increase of the current suppresses this feature, leading to a more simple picture of inverted 
MIS oscillations. 

III. THEORY 

The theoretical interpretation of our data is based on the physical model of Dmitriev et 
al.J generalized to the two-subband case. The elastic scattering of electrons is assumed to 
be much stronger than the inelastic one. This scattering maintains nearly isotropic carrier 
distribution at moderate currents, when the momentum gained by an electron moving in 
the electric field between the scattering events is much smaller than the Fermi momentum. 
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FIG. 3: (Color online) Magnetoresistance of the sample B at T = 1.4 K. The values of the current 
are 10 (bold), 30, 50 (dash), 100 (bold dash), 150, 200 (short dash), and 300 (bold) /uA. The inset 
shows amplitudes of the inverted peaks at B = 0.34 T. 

Since the intersubband elastic scattering is also much stronger than the inelastic scattering, 
the isotropic part of electron distribution function, Z^, is common for both subbands and 
depends only on the electron energy e. When the current of density j flows through the 
sample, the kinetic equation for this function is written as 



P d d 
D^a^ oe oe 



(1) 



where P = j'^pd is the power of Joule heating (the energy absorbed per unit time over a 
unit square of electron system) expressed through the diagonal resistivity p^, and is the 
density of states. The function (yd{s) can be written through the electron Green's functions, 
which are determined by the interaction of electrons with static disorder potential in the 
presence of magnetic field. The free-electron states in the magnetic field described by the 
vector potential (0, -Bx, 0) are characterized by the quantum numbers j, n, and Py, where 
j = 1,2 numbers the electron subband of the quantum well, n is the Landau level number, 
and Py is the continuous momentum. Using the free-electron basis, one obtains 
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FIG. 4: (Color online) Magnetoresistance of the sample B at T = 4.2 K. The values of the current 
(/uA) are 1, 50, 80, 100, 120, and 150 for the curves marked by the numbers from 1 to 6, respectively. 
The other curves corresponds to the currents of 200 (short dash), 250 (bold dash) 300 (solid), 350 
(dash), and 400 (bold) ^A. 

where e is the electron charge, m is the effective mass of electron, G^-^''** are the retarded 
(s = R) and advanced (s = A) Green's functions, and is the normalization square. The 
Zeeman splitting is neglected, so the electrons are assumed to be spin-degenerate. The 
double angular brackets in Eq. (3) denote averaging over the random potential. In terms of 
the Green's functions, the density of states is given by 

jnpy j 

The dimensionless function Sj^ is found from the implicit equation 

271 ^ e- hwcin + 1/2) - Sj - Sj^' ^ ' 

J, Jf 

where Uc is the cyclotron energy, ej is the subband energy, and Tjji are the quantum lifetimes 
of electrons with respect to intrasubband {j' = j) and intersubband {j' ^ j) scattering. 
Equation (5) is valid when the correlation length of the disorder potential is smaller than 
the magnetic length, and the disorder-induced energy broadening of the subbands is smaller 
than the subband separation A12 = 62 — si. It corresponds to the the self-consistent Born 
approximation (SCBA). 



According to the definition (2), the diagonal conductivity is 

(7d = y de (^d{£)- (6) 

Therefore, multiplying the kinetic equation (1) by the density of states D,, and energy e, and 
integrating it over £, one obtains the balance equation P — Pph, where Pph — — J deeD^Js{f) 
is the power lost to the lattice vibrations (phonons). 

Below we consider the case of classically strong magnetic field, uJcTtr ^ 1, when (Td{s) is 



written in terms of Sj^ as 
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where ni and n2 are the electron densities in the subbands, Ug — ni + n2, and rjj, are the 
transport times of electrons. Both Tjf and Tjj, are determined by the expressions 

where Wjji{q) are the Fourier transforms of the correlators of the scattering potential, 
Fjj'{9) = 1 — 2kjkjf cos 9/ {k'j + and kj is the Fermi wavenumber for the subband j. 
The electron densities in the subbands are expressed as rij — kj/2'!r. 

In DQWs, where the energy separation between the subbands is usually small com- 
pared to the Fermi energy, the difference kf — /c^ is small in comparison with kf + /cf so 
that rii ~ 77,2 — ns/2. Furthermore, in the symmetric (balanced) DQWs, where the elec- 
tron wave functions are delocalized over the layers and represent themselves symmetric and 
antisymmetric combinations of single-layer orbitals, one has nearly equal probabilities for 
intrasubband and intersubband scattering owing to Wu{q) ~ ^22(5) — w^iq), provided that 
interlayer correlation of the scattering potentials is weak. Therefore, Tjj ~ T12 ^ 2t, and 
Tjj ~ — 2Ttr, where r and Ttr are the averaged quantum lifetime and transport time, 
respectively. In these approximations, Eq. (7) is written in the most simple way: 

ad{e)c^afvl = ^(2^ie + 2^2e), Vj, ^ 21mSj, (9) 

where a'^^ — aj^po, cr± — e^ns/mujc is the Hall conductivity, and po = m/e^Ttrns is the 
classical resistivity. The function T>^ — is the dimensionless density of states, containing 
oscillating (periodic in tvjJc) part 7^. Therefore, it is convenient to solve the kinetic equation 
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by representing the distribution function as a sum /° + Sf^, where the first term slowly 
varies on the scale of cyclotron energy, while the second one rapidly oscillates/ The first 
term satisfies the equation 

(10) 

Solution of this equation can be satisfactory approximated by a heated Fermi distribution. 
This is always true if the electron-electron scattering dominates over the electron-phonon 
scattering and over the electric-field effect described by the left-hand side of Eq. (10). In 
this case, the Fermi distribution of electrons is maintained against the field-induced difi^usion 
in the energy space, while the electron-phonon scattering determines the effective electron 
temperature Tg. In the general numerical solution of Eq. (10) involving electron- 

phonon scattering in the coUision integraP^ confirms that /° is very close to the heated 
Fermi distribution. 

The equation for the oscillating part, 5fs, is then written in the following form: 

Below we search for the function 5fs in the form Sfe — {df^ /de)(fe, where is a periodic 
function of energy. Taking into account that the main mechanism of relaxation of the 
distribution 5fs is the electron-electron scattering, one may represent the linearized collision 
integral Jei^Sf) as 

1 df^ 1 

X[ip, + ip,, -ip,+5,-ip,,_Se])e',5e^ ^ ^ ^33',h3'v (12) 

where Se is the energy transferred in electron-electron collisions, Mjj'j^j/^ is the probability 
of scattering (when electrons from the states j and j' come to the states ji and ][)■, M is the 
normalization constant, and the angular brackets (. . .)^, denote averaging over the energies 
e' and 5e. Expression (12) is a straightforward generalization of the result of Ref. 7. The 
characteristic inelastic scattering time Ti„ describes the relaxation at low magnetic fields, 
when Vj^ are close to unity. In this case the collision integral acquires the most simple form 
Je{Sf) = —Sfe/Tin, i.e. the relaxation time approximation is justified. 
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The resistivity pa — c^^V^^l written, according to Eq. (6), in the form 

.. = ,/...;(-f)(i.|.), (13) 

where we have taken into account that dfe/de ~ {df'^/de)[l + dipe/de\. Therefore, in 
order to calculate the resistivity, one should find iff, by using Eqs. (11) and (12). In 
general, Eq. (12) is an integro-differential equation that cannot be solved analytically. 
However, the property of periodicity allows one to expand in series of harmonics, = 
exp(27ri/c£//ki;c), and represent Eq. (11) as a system of linear equations: 

oo t 

(g-^ + A;2)(^fc+ [{2kk'-k'^)^k-k' + Q-'Ckk']'Pk'^'2ik^^k: (14) 

fc'=— oo 

where 

Q-^J- (15) 

is a dimensionless parameter characterizing the nonlinear effect of the current on the trans- 
port. The matrix Ckk', whose explicit form is not shown here, describes the effects of 
electron-electron scattering beyond the relaxation time approximation. 

The harmonics of the density of states, 7^, as well as the coefficients Ckk', which are 
expressed in terms of products of these harmonics, are proportional to the Dingle factors 
exp(— /ctt/cUct). Therefore, searching for the coefficients ipk at weak enough magnetic fields, 
when e~'^/^'='^ is small, one can take into account only a single {k = ±1) harmonic. Within 
this accuracy, one should also neglect the sum in Eq. (14). This leads to a simple solution 
(fi±i = ±ij±i{hujjTr)Q/{l + Q). Since 7+1 7_i = -2e-''/'^='^ cos(7rAi2//ia;c), Eq. (13) is 
reduced to a simple analytical expression for the resistivity: 

^ = l + e--/-i^ fl + cos^ 

The second term in this expression, proportional to g^^'^/'^'^^, differs from a similar term 
of the single-subband theory^ by the modulation factor [1 + cos(27rAi2//jti;c)] /2 describing 
the MIS oscillations. The last term in Eq. (16) describes the SdHO, which are thermally 
suppressed because of the factor T = (27r^Te/ hujc) / sinh(27r^Te/ hwc) ■ The Fermi energy is 
counted from the middle point between the subbands, {ei + e2)/2, and, therefore, is directly 
proportional to the total electron density, Sp — h^Tms/^m. 
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IV. RESULTS AND DISCUSSION 



The basic features of our experimental findings can be understood within Eqs. (16) 
and (15). In the hnear regime, when the parameter Q is small, this equation gives a good 
description of the MIS oscillations experimentally investigated in Ref. 10. As the current 
increases, the amplitudes of these oscillations decrease, and then the flip occurs, when the 
MIS peaks become inverted. In contrast, the SdHO peaks are not affected by the the current 
directly, and their decrease is caused by the effect of heating. The flip of the MIS oscillations 
corresponds to Q — 1/3. Since Q is inversely proportional to the square of the magnetic 
field, there exists the inversion field, determined from the equation Q — 1/3, where Q 
is given by Eq. (15). This feature is observed in our experiment, see the inset to Fig. 2. For 
the sample B, we have extracted Binv for several values of the current. The results are shown 
in Fig. 5. At 4.2 K the experimental points follow the linear Biny{I) dependence predicted 
by Eq. (15). Since the ratio Bin^/I is proportional to the square root of the inelastic 
relaxation time Tin, we are able to estimate this time from experimental data as Tin — 64 
ps at T = 4.2 K. Assuming the T"^ scahng of this time,^ one obtains h/Tin — 6.6 mK at 
T = 1 K, which is not far than the theoretical estimate h/Tin — 4 mK at T = 1 K based 
on the consideration of electron-electron scattering.''' The positions of experimental points 
at T = 1.4 K also fit within this picture if the electron heating is taken into account. The 
increase of electron temperature with increasing current (heating effect) leads to deviation 
of the Binv{I) dependence from linearity because of temperature dependence of Tin, and this 
deviation is essential at T = 1.4 K; see Fig. 5. The same consideration, applied to the 
high-mobihty sample A, gives the inelastic scattering time Tj„ ~ 108 ps at T = 4.2 K, which 
is very close to the theoretical estimate. 

When the current becomes high enough (Q ^ 1), Eq. (16) predicts saturation of the 
resistance, when the amplitudes of inverted MIS peaks are three times larger than the 
amplitudes of the MIS peaks in the linear regime {Q <^ 1). We indeed observe the regime 
resembling a saturation, with almost three times increase in the amplitudes of inverted peaks 
for both samples at T = 1.4 K (see Figs. 1 and 3). For higher temperatures the behavior is 
similar, though the maximum amplitudes of inverted peaks are only slightly larger than the 
amplitudes in the linear regime. We explain this by the effect of heating on the characteristic 
times. Though the resistivity in the high-current regime {Q ^ 1) no longer depends on Tin, 
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FIG. 5: (Color online) Dependence of the inversion field on the current for the sample B at T = 4.2 
K and T = 1.4 K. (points) The dashed lines correspond to a linear Binv{I) dependence assuming 
Tin = 64 ps at 4.2 K (580 ps at 1.4 K). The solid lines represent the calculated Binv{I) dependence 
taking into account electron heating by the current. 

there is a sizeable decrease in the quantum lifetime r with increasing temperature/^ which 
takes place because the electron-electron scattering contributes into r. As a result, the 
Dingle factor decreases, and the quantum contribution to the resistance becomes smaller 
as the electrons are heated. At higher initial temperature, when Tin is smaller, the regime 
Q ^ 1 requires higher currents. The corresponding increase in heating reduces the quantum 
contribution, so the maximum amplitudes of inverted peaks never reach the theoretical limit 
and are expected to decrease with increasing initial temperature. The slow suppression of the 
inverted peaks with further increase in the current (see the inset to Fig. 3) is explained by 
the same mechanism. This conclusion is supported by the experimental observation that the 
suppression is more efficient at lower magnetic fields, when the Dingle factor exp(— tt/cUct) 
is more sensitive to the temperature dependence of quantum lifetime r. 

To illustrate the above- discussed relation of the basic theoretical predictions to our ex- 
periment, we present the results of theoretical calculations according to Eqs. (15) and (16) 
in Fig. 6. The calculations are done for the sample B at 4.2 K, so the theoretical curves 
show the expected behavior of the measured magnetoresistance from Fig. 4. We take into 
account the effect of heating, described by using the collision integral for interaction of 
electrons with acoustic phonons^^ and temperature dependence of the quantum lifetime r 
of electrons determined empirically from the studies of the MIS oscillations in the linear 
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FIG. 6: (Color online) Calculated magnetoresistance of the sample B at T = 4.2 K and different 
currents: 1, 50, 80, 100, 120, and 150 /^A for the curves marked by the numbers from 1 to 6; the 
other curves corresponds to I = 200 (short dash), 250 (bold dash), 300 (solid), and 400 (bold) 
/iA. The additional (dashed) line 1 shows the linear magnetoresistance determined by the SCBA 
calculation of the density of states in Eq. (13). 

regime. The theoretical plots demonstrate a reasonable qualitative agreement with the ex- 
periment. However, the theory predicts a slower suppression of the inverted peak amplitudes 
with increasing current at weak magnetic fields. This may be a consequence of underesti- 
mated heating, because the screening effect on the electro n-phonon interaction^'^ has not 
been taken into account in the calculation of the power loss to acoustic phonons. Similar 
calculations carried out for different samples at different temperatures are also in agreement 
with experimental data. 

The simple theory fails do describe the interesting and unexpected feature observed in 
our experiment, the current-induced splitting of the MIS oscillation peaks. This kind of 
nonlinear behavior is well-reproducible, we see it in different samples. We have found that 
a possible explanation of this feature can be based on the theory presented in Sec. Ill, if 
higher harmonics of the distribution function Sf^ are taken into account. We have carried 
out a numerical solution of the system of equations (14) under some simplifying assumptions 
about the collision integral. In the first case, we have assumed equal probabilities for all 
possible electron-electron scattering processes, so the matrix Mjji^^j'^ in Eq. (12) is replaced 
by a constant. Another limiting case we consider is the complete neglect of intersubband 
transitions in electron-electron collisions, when Mjj'j-^j/^ oc ^jjx^j'j^- This case is also rea- 



13 



0.8 I 1 ' 1 ' 1 ' 

0.2 0.4 0.6 0.8 
1.4, 




0.8 I ' ' ' ' ' ' ' ^ 

0.2 0.4 0.6 0.8 
B(T) 

FIG. 7: (Color online) (a) Calculated magnetoresistance of the sample B at T = 4.2 K and I = 120 
fiA. The plot 1 correspond to simple theory [Eq. (16)], while the others represent the results of 
numerical solution of Eq. (14) for the cases of subband-independent electron-electron scattering 
(2) and only intrasubband electron-electron scattering (3). (b) The same plots, where the SdHO 
contribution is excluded. The density of states is found within the SCBA. 

sonable, since electron-electron scattering at low temperatures assumes a small momentum 
transfer, so the intersubband scattering contribution should be suppressed owing to reduc- 
tion of the overlap integrals of envelope wave functions of electrons. Then, the coefficients 
and Ckk' have been determined by using the density of states numerically calculated within 
the SCBA; see Eq. (5). The results, corresponding to / = 120 fiA for the sample B are pre- 
sented in Fig. 7. In the low-field region, where the MIS peaks are inverted, the calculation 
shows a considerable increase in their amplitudes above 0.2 T, where contribution of higher 
harmonics of the density of states becomes essential. This enhancement occurs because of 
the current-induced mixing between different harmonics of the distribution function, for- 
mally coming from the term with '~fk-k' in the sum in Eq. (14). In contrast, in the linear 
regime, the SCBA magnetoresistance is close to the magnetoresistance calculated within 
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the single-harmonic approximation [Eq. (16)]; see Fig. 6. Above 0.27 T, where the Landau 
levels become separated, one can see features associated with the specific semi-elliptic shape 
of the SCBA density of states. In the vicinity of the inversion field (Bi^v — 0.4 T), where 
the contribution of the first harmonic of the distribution function is suppressed {Q ~ 1/3) 
while the higher harmonics are still active, two sets of MIS peaks arc seen. It is not surpris- 
ing, because higher harmonics of the density of states contain the factors cos(/i;7rAi2//k<;c) 
describing higher harmonics of the MIS oscillations. Above the inversion field, the resistance 
is considerably smaller than the resistance predicted by the single-harmonic approximation, 
and a sphtting of the MIS peaks occurs. The splitting increases with the increase of the 
magnetic field. These effects are caused by the contribution of higher harmonics of the 
density of states in the collision integral. Indeed, in the single-harmonic approximation the 
collision integral contains only the outcoming term proportional to ip^. This approximation 
becomes insufficient in higher magnetic fields, when incoming terms in the collision integral 
(12) are also important, so the relaxation of the distribution function, which counteracts 
the diffusion of electrons in the energy space, becomes less efficient. This means that the 
effect of the current on the distribution function increases, and the resistance is lowered. 
The described suppression of the collision-integral term is more significant in the regions of 
the MIS resonances, when Au/hu is integer, because the peaks of the density of states are 
the narrowest in these conditions, and the energies transferred in the electron-electron col- 
lisions, Se, are small. Away from the MIS resonances, the energy space for electron-electron 
scattering increases, especially when the intersubband transitions are allowed (see curve 2 
in Fig. 7). Therefore, the relaxation is less suppressed as compared to the center of the 
MIS peak, and the effect of the current is weaker. The above consideration explains why 
the centers of the MIS peaks drop down, so the peak splitting takes place. 

The SCBA has a limited applicability for description of the density of electron states in 
the magnetic field. In particular, it leads to non-physically sharp edges of the density of 
states, which generate the harmonics with large k in Eq. (14). This apparently leads to 
an overestimate of the effect of the current on the resistance in the region where the MIS 
peaks are inverted, see Fig. 7. To avoid such singularities, and to have a further insight into 
the problem of nonlinear magnetoresistance, we have considered the expression 
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FIG. 8: (Color online) The same as in Fig. 7 for the Gaussian model of the density of states. 

which corresponds to the Gaussian model for the density of states and describes two inde- 
pendent sets of Landau-level peaks from each subband (strictly speaking, the Landau-level 
peaks are not independent because of elastic intersubband scattering, as follows from Eq. 
(5), see more details in Ref. 14). The magnetic-field dependence of the broadening en- 
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ergy F has been set to make the first [proportional to cos(27r£:/^c)] harmonics of T) 



and equal. The results of the calculations using V^-^ instead of the SCBA density of 
states are shown in Fig. 8. The magnetoresistance in the region of inversion appears to 
be nearly the same as predicted by the simple single-harmonic theory. In the region above 
the inversion field, the splitting of the MIS peaks does not take place if the intersubband 
electron-electron scattering is forbidden. This is understandable from the discussion given 
above: if different subbands contribute into the density of states independently, the effi- 
ciency of electron-electron collisions does not depend on the ratio /hjj and the reduction 
of the collision integral owing to incoming terms causes just a uniform suppression of the 
whole MIS peak. In the SCBA, when the shape of depends on this ratio, the splitting 
of the MIS peaks does not necessarily require the intersubband electron-electron scattering. 
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FIG. 9: (Color online) Evolution of the nonlinear magnetoresistance calculated using the param- 
eters of the sample B when the current varies from 100 to 150 /iA with the step of 10 /lA. The 
Gaussian model of the density of states and the assumption of subband-independent electron- 
electron scattering are used. 

If the intersubband electron-electron scattering is allowed, the magnetoresistance pictures 
obtained within the Gaussian model, as well as within the SCBA model above the inver- 
sion point, qualitatively reproduce the features we observe experimentally. The results of 
calculations presented in Fig. 9 demonstrate that varying the current in a relatively narrow 
range leads to a dramatic reconstruction of the magnetoresistance oscillation pattern. 




FIG. 10: (Color online) Comparison of the measured and calculated nonlinear magnetoresistance 
in the sample A at T = 4.2 K and I = 75 /xA. The Gaussian model of the density of states and the 
assumption of subband-independent electron-electron scattering are used in the calculations. 



Numerical calculation of magnetoresistance in the high-mobility sample A also gives the 
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results very similar to what we see experimentally. To demonstrate this, we have put exper- 
imental and calculated curves together in Fig. 10. Apart from a weak negative magnetore- 
sistance at low fields and a slight decrease in the MIS oscillations frequency with increasing 

B (the features we see in all our samples^^'^^ both in linear and nonlinear regimes), the 
agreement between experiment and theory is good. 

V. CONCLUSIONS 

Investigation of nonlinear transport of 2D electrons in magnetic fields enriches the knowl- 
edge of the quantum kinetic properties of electron systems and of the microscopic processes 
responsible for the observed modifications of the resistivity. In our work, we have demon- 
strated that using double quantum well systems opens wide possibilities for studying the 
nonlinear behavior. The presence of the MIS oscillations, which modulate the quantum 
component of the resistivity, allows us to investigate the current dependence of the quan- 
tum magnetoresistance. In particular, we are able to determine the magnetic fields Bin^ 
corresponding to the current-induced inversion of the magnetoresistance. This inversion 
manifests itself in a spectacular way, as a fiip of the MIS oscillation pattern. We point out 
that this behavior resembles recently obscrvcd^^ inversion of the MIS oscillations by the 
low-frequency (35 GHz) microwave radiation. This is not surprising, because the physical 
mechanism in both cases is similar. Apart from the fiip of the MIS oscillations, we have 
observed a wholly unexpected quantum phenomenon, the splitting of the MIS oscillation 
peaks in the region of fields above the inversion point 

We have shown that the theoretical explanation of all the observed phenomena can be 
based on the kinetic equation for the isotropic non-equilibrium part of electron distribution 
function. This function oscillates with energy owing to oscillations of the density of electron 
states in the magnetic field. The effect of electric current on this function, the increase 
of electron diffusion in the energy space, is equilibrated by the inelastic electron-electron 
scattering. Theoretical explanation of the most of observed phenomena is done in a simple 
single-harmonic approach, which allowed us to determine the inelastic relaxation time Tin 
by comparison of experimental data with theory. The values of rj„ for different samples 
are close to the theoretical estimates of this time, and confirm the predicted^ temperature 
dependence Tin T~^. Thus, our data on the inelastic relaxation time in double quantum 
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well samples are in agreement with the data obtained in single quantum well samples.^ 
The description of the sphtting of MIS oscillations requires a more detailed numerical 
analysis including consideration of higher harmonics of both the density of states and the 
distribution function. Apart from the verification of the basic principles of the theory of 
Ref. 7, this analysis demonstrates sensitivity of the nonlinear behavior to the shape of the 
density of electron states and to the details in description of inelastic scattering. Therefore, 
investigation of nonlinear magnetoresistance in relatively weak magnetic fields offers a 
tool for studying the electron states and scattering mechanisms both in single and double 
quantum wells. 
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